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ABSTRACT 

Motivated by the dawn of preeision eosmology and the weahh of forthcoming high pre- 
cision and volume galaxy surveys, in this paper we study the effects of inhomogeneities 
on light propagation i n a flat ACDM backgro und. To this end we use exact solutions of 
Einstein's equations (iMeures fc Brunill2011l ) where, starting from small fluctuations, 
inhomogeneities arise from a standard growing mode and become non-linear. While 
the matter distribution in these models is necessarily idealised, there is still enough 
freedom to assume an arbitrary initial density proflle along the line of sight. We can 
therefore model over-densities and voids of various sizes and distributions, e.g. sin- 
gle harmonic sinusoidal modes, coupled modes, and more general distributions in a 
ACDM background. Our models allow for an exact treatment of the light propaga- 
tion problem, so that the results are unaffected by approximations and unambiguous. 
Along lines of sight with density inhomogeneities which average out on scales less than 
the Hubble radius, we flnd the distance redshift relation to diverge negligibly from the 
Friedmann-Lemaitre- Robertson- Walker (FLRW) result. On the contrary, if we observe 
along lines of sight which do not have the same average density as the background, 
we find large deviations from the FLRW distance redshift relation. Hence, a possibly 
large systematic might be introduced into the analysis of cosmological observations, 
e.g. supernovae, if we observe along lines of sight which arc typically more or less dense 
than the average density of the Universe. In turn, this could lead to wrong parameter 
estimation: even if the Cosmological Principle is valid, the identification of the true 
FLRW background in an inhomogeneous universe maybe more difficult than usually 
assumed. 

Key words: cosmology: theory, gravitation, large-scale structure of Universe, dark 
energy, methods:analytical, supernovae: general 



1 INTRODUCTION 

At the very basis of modern cosmology lies the assump- 
tion the Universe is, at any given time, homogeneous and 
isotropic on large scales. This is translated mathematically 
into a Robertson- Walker metric, i.e. a metric that is assumed 
to represent a space average and is therefore exactly homoge- 
neous and isotropic. In addition, the non-trivial hypothesis 
is made that this metric should be a solution of Einstein's 
equations, thereby giving rise to a FLRW universe model. 
Having assumed the Cosmological Principle, the growth of 
inhomogeneities and their effects are typically modelled with 
perturbation theory about a "background" FLRW model. 
Within this framework, the formation of non-linear struc- 
tures at smaller scales is considered in N-body simulations 
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using the Newtonian approximation IjSpringel et al.ll2005l ). 
Most observations are interpreted assuming this Friedman- 
nian framework; in particular, distances are computed as- 
suming a FLRW distance-redshift relation, i.e. completely 
neglecting inhomogeneities. 

We are currently living in a time when galaxy surveys 
and other observations are reaching unprecedented sky cov- 
erage and precision; therefore, it seems timely to fully in- 
vestigate the effects of the non-linear growth of structures 
on observations, within a general relativistic framework. 
Much work has been done in trying to understand the effect 
of inhomogeneities on ob servations by us ing perturb ations 
around an FLRW model llDver fc Roedeij[l9 72. 1973 1 Il974 
ISas akil[l987l: iFutamase fc Sasakilll989l: IPvne fc Birkinshav 



2004; iBarausseTt al.l l2005l : iBonvin et all l2006l : iRasanen 
120111 ). However, these analyses are limited to linear 
structure growth and therefore cannot properly take 
into account non-linear inhomogeneities. Non-linearities 
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can be modelled using idealised matter distributions. 
Several different approaches using generalised Swiss- 
Cheese models or alternative geo metrie s have been 



consid ered 



m 



e.g. iBrouzakis et al] toOTl): iMarra et all 



Brouzakis et al.l (|2008|'1: iBiswas fc Not aril ll2008h 



IChfton fc Zuntzl (|2 0091 : Bo leiko fc CelerieJ (|201ol ): ISzvbkal 



l 2011f ): lNwankwo et al. (20f 1 ). where significant deviations 



are usually only found for very large scale density inhomo- 
geneities. On the contrary, claims are made that small scale 
structure formation might have a back- reaction e f fect o n the 
overall expansi on of the Universe, see iRasanenI (|2006l ) and 
iBucherd l|2008h for overviews of this topic. We do not in- 
vestigate this phenomenon of back-reaction in this paper 
though, as our model clearly splits into inhomogeneities and 
background dynamics. 

A class of models which are very appropriate for con- 
sidering a discrete distribution of matter in an otherwise 
FLRW e xpanding universe has re cently been analysed in 
depth bv lClifton fc Ferr eira (2009a b) . In these mo dels, orig- 
inally in troduced b y Lindquist fc Wheelen (|l957l ) and re- 



vised bv iRedmountI l| 19881 ). matter is described by pointlike 
masses in a spherically symmetric void box (represented by 
Schwazschild space-times) and these boxes are distributed in 
a lattice and the overall expansion is described by the Fried- 
mann equation. The main motivation of Clifton and Ferreira 
in following the Lindquist and Wheeler construction is the 
observation that the Universe largely consists of galaxies 
and clusters of galaxies surrounded by vacuum. The ques- 
tion they address is how observations and measurements of 
the cosmological parameters are affected in a highly inhomo- 
geneous universe whose overall dynamics are homogeneous 
and isotropic. However, this lattice construction is only an 
approximate solution to Einstein's equations and has regions 
of "no man's land" in between the matche d sphere s whic h 
might have an effect on the light tracing, see lCliftonI (|201l|). 
Perhaps even more importantly, the inhomogeneities in this 
model are strongly non-linear at all times. 

Having the s ame t ype of questions addressed by 
IClifton fc Ferreiral l|2009al ) in mind, in this paper we shall 
investigate the optical properties of an exact solution to 
Einstei n's field equations (EF Es), developed in our previous 
paper (Mcuros & Bruni 201l]) . This solution belongs to the 
second clas s Szckcrcs models (|Szekereslll975l ) for dust, gen- 
eralised by iBarrow fc Stein-Schabes (119841) to inc lude the 
cos mological constant A. InlMeures fc Brunil (|201ll ) we used 
the ICoode fc Wainwrightl (|l982al ) formulation of the Szek- 
eres models to reconsider the Barrow and Stein-Schabes so- 
lutions in a ACDM context. In our model, starting from 
standard small perturbations of a ACDM homogeneous and 
isotropic universe, the matter distribution is continuous and 
can evolve to a highly non-linear stage. In the process, the in- 
homogeneities can either form a distribution of large voids, 
over-densities, or a mixture of the two with over-densities 
possibly even forming pancakes as in the Zel'dovich approx- 
imation in Newtonian cosmology. The benefit of our model 
is therefore two- fold: i) we consider exact solutions of Ein- 
stein's equations, therefore avoiding any possible problem 
associated with approximations and matching and ii) these 
exact solutions describe non-linear inhomogeneities grow- 
ing on top of a FLRW background with the possibility of 
modelling a rather arbitrary distribution of both voids and 
over-densities. 



In particular, we shall consider the effect of inhomo- 
geneities on the redshift, angular diameter distance and 
distance modulus. We begin by considering single mode 
harmonic sinusoidal deviations from homogeneity and then 
the case of coupled r nodes. Finally, closer in spirit to 
the IClifton fc Ferreiral (|2009aj) work, we consider inhomo- 
geneities where peaks and voids, arranged in a periodic ar- 
ray, are more strongly separated than can be achieved by 
simple harmonic distributions. We demonstrate that the de- 
viations from the FLRW background in the determination 
of the distances is mainly due to the Ricci and Weyl fo- 
cusing terms in the Sachs equations and show that instead 
the shear of the null congruence has a negligible effect. We 
also briefiy investigate the effect of mode coupling on the 
growth of structure and, interestingly, we show that even 
a long wavelength mode with small amplitude can strongly 
enhance the growth of short wavelength modes, thanks to 
the non-linear coupling. The non-linear interaction of differ- 
ent modes does not seem to influence the distance measures 
significantly and we find that changes in the redshift and 
distances are mostly affected by the long wavelength modes. 
Using an array of density profiles which are not sinusoidal 
but quite peaked around the maximum and separated by 
large voids, we find that the effect on the redshift and dis- 
tance measures does not prove to be significantly more than 
using an initially sinusoidal density distribution with the 
same wavelength of the array scale. Overall, we find that 
all deviations in the redshift and distance measures are less 
than 1%, when we consider what we refer to as "compen- 
sated inhomogeneities" along the line of sight, i.e. where the 
average density along the line of sight matches the back- 
ground density. However, this does not need to be the case: 
when the inhomogeneities are on average above or below the 
background, the effects on redshift and distance measures 
can be very large. 

A summary of the paper is as follows. In Sec. [51 we 
will briefly present the exact solu tion we will be u s ing in 
this work, referring the reader to iMeures fc Brunil (|201l|) 
for more details. Subsequently, in Sec. (31 we will derive the 
null geodesic equations for the given metric and in Sec. |4]we 
derive the form that the Sachs optical equation will take for 
a given physical situation that is to be investigated. In Sec. 
[5] we shall present the results of our analysis, considering 
single mode deviations in Sec. 15. II multiple modes and their 
coupling and effects in Sec. 15.21 and the effects of an array of 
strong peaks and large voids in Sec. 15.31 In Sec. [6] we draw 
our conclusions. In Appendix |X| we present details on the 
tetrad transformations needed to derive the Sachs optical 
equations in our model. 

Throughout the paper we choose units c = SnG = 1 
and assume the standard JIa = 0.75 and in the commonly 
used units Ho — 72kms~^Mpc~^. 



2 AN EXACT SOLUTION OF SUFFICIENT 
GENERALITY 

In this section, we would like to present a sho rt summary 
of th e main results of our previous paper (Me ures fc Brunil 
I2OIH ) , where we developed a class of exact solutions to EFEs, 



Gab ~ Tab — ^Qab, (1) 
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for the metric 



ds 



'dt^ + S{tf [dx'^ + dy"^ + Z(x, tfdz^] , (2) 



where we are using synchronous comoving coordinates, so 
that t represents a universal cosmic time. The variable x 
denotes all three spatial Cartesian coordinates x, y and z. 
Note that whenever we will use the variable z, we will be 
referring to the coordinate; whereas we will refer to the red- 



shift as 



respectively for an inhomogeneous 



or FLRW universe. We solve Eq. ([T} for an irrotational pres- 
sureless fluid with four velocity u"" = Sq, for which 



T = p, 



(3) 



and all other components of T*^" are zero. The so- 
lutions for a metri c of this type for pure dust were 
first introduced by ISzekered (119751 ) and then brought 
into a simil ar notation that we are us ing here and 
analysed by ICoode fc Wainwrighd (jl982al lbh. The solu- 
tion for dust and a cosmologi cal constant was found by 
iBarrow fc Stein-Schabej ([198J). We are here considering 
one specific sub-class of these models, usually called the 
second class Szekeres models, and we choose the FLRW 
background to be spatially flat. For de tailed accounts o f 
exact solutions in ge neral relativity, see Krasinska ([1993), 
IStephani etHI (|2003l ) and lBoleiko et all (|2009l ). 

The function S{t) is the scale factor of a FLRW ACDM 
background, with Friedmann equation 



which admits the solution 



353 3 ' 



S{t) = 



!^A 



1/3 



sinh 



2/3 



-HoV^t 



(4) 



(5) 



Here p — po/S^ is the background energy-density and we 
have normalised S — 1 today. We have used the standard 
parametrisation Q,m ~ po/(3-ffo) and JIa = A/{3Hq), where 
Ho is the Hubble parameter and Sim = 1 — ^a- 

Note that we choose S as well as the function ^ in (l2|) to 
be dimensionless, so that the coordinates have a dimension 
of length. 

The function Z in the line element ^ can be split as 

Z{x,t)^F{z,t)+A{x), (6) 

where A can be written in the form 

yl(x) = 1 + Bp+{z) {[x + 7(2)]^ + [y + u^{z)f} , (7) 
where 

(8) 



B^^HilnAii-nAff'., 



with the free functions l3+(z), 7(2) and lj{z). 

Remarkably, F obeys the second order linear homoge- 
neous ordinary differential equation 



F + 2^F-^F = 0, 



(9) 



which is exactly the same equation that 5 = {p— p)/p satis- 
fies in Newtonian linear perturbation theory, as well as in the 
synchronous comoving gauge for relativistic perturbations of 




2 3 4 5 




Figure 1. Plots of the growing (top panel) and decaying (bottom 
panel) modes of the solution for F, as derived from Eq. l[9}. The 
solutions plotted here are given in Eqs. mial l and l lllbl l. 



FLRW. This equation admits two linearly independent so- 
lutions and hence we write 



F{z,t)^l3+{z)f+{t) + l3.{z)f-{t), 



(10) 



where /3+ and /3_ are free functions of the coordinate z and 
/+ is the growing mode and /_ is the decaying mode of the 
solution, which we find to be 



/- 
/+ 



cosh(r) 
sinh(r) 
cosh(r) 



sinh2/3( 



sinh(r) J cosh^(r^ 



-dr, 



(11a) 
(lib) 



where we defined the dimensionless variable r = |_ffoV"Ai. 
The two independent solutions for F are shown in Fig. ^. In 
the matter dominated era, when the effect of A is negligible, 
f+j t) PC S(t ) PC t^ '^ and /_ ex i~^, as it is well known, see 
e.g. |Peeblesl(|l980l ). 

With all the free functions having been identified and 
having solved for all the time dependent functions, we can 
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find the expressions for the density p, the background den- 
sity p and density deviation 5 to be 



P = 



and 



5 = 



poA 



S3{F + A)' 
53' 



P 
F 
Z 



F 



F + A' 



(12) 
(13) 

(14) 
(15) 



As we are only concerned with late times here, we will not 
be considering the decaying mode of F, /3_/_, as it would 
not effect our results, but complicate the calculations. This 
can easily be done by choosing /3_ = in (|10|l . On the 
other hand, a very useful feature of our solution is that F 
satisfies the linear differential equation ((9)l, so that a super- 
position principle applies, and that the arbitrariness of the 
growing mode function /?+ in (|1U|I allows us to construct an 
arbitrary matter distribution along z. We find that at early 
times, along the z-axis, S ~ —F and so choosing a func- 
tion /3+ directly determines the initial matter distribution 
along the z-axis. Choosing /3+ = Asm{kz) and therefore 
S oc sin(fc2:) initially for some amplitude A and some wave- 
length k, implies that /3+ and hence the initial distribution 
of 5 are periodic on a scale of 2ii /k. This situation is what 
we will be referring to as a compensated density deviation 
later on in the paper, as averaging the density deviation 5 
along the «-axis at early times, would tend to a zero average 
once the period of the deviations is reached. In comparison, 
we will be referring to over-densities for /3+ < 0, i.e. 5 > 
and under-densities for /?+ > and hence 5 <Q. 

To aid the reader in gaining an intuitive understanding 
of the density profiles that we will be using in this paper, 
we have included two figures with density profiles for 7 = 
a; = 0. In Fig. [2] we show what the shape of the density 
deviation would be today, if we chose the initial profile to 
be (5 oc sin(fcz;) for k = 27r/8Mpc~^ along the z-axis. Only 
two spatial dimensions are displayed here, because of the 
symmetry of the model around the z-axis for 7 = a; = 0. 
Note that in Figs. [2l [3l and [TOl we use coordinate distances, 
measured in Mpc, which serve as a reference distance. Given 
our metric ([2]) and our choice S = 1 today, the coordinate 
distance in directions orthogonal to the z-axis equals the 
proper distance today, but in any other direction the proper 
distance will involve an integration over the function Z. 

In Fig. Owe display how an initial sinusoidal density de- 
viation grows non-linearly into a shape with voids and high 
over-density peaks. While the metric function F, initially 
F oc S, evolves linearly and remains sinus oidal, 5 can grow 
strong ly, even developing pancakes, see iMeures fc Brunil 
(|201ll '). 



3 THE NULL GEODESIC EQUATIONS 

We now want to calculate the null geodesic equations us- 
ing the Euler-Lagrange formalism. We start from the La- 




r [Mpc] 



z [Mpc] 



Figure 2. Density deviation profile today, corresponding to an 
initial density perturbation of <5 oc sin(fc2) for k = 27r/8Mpc~^ 
along the z-axis, for 7 = oj = 0. We only indicate the coordinate 



distance r = \J x^ -\- y^ from the z-axis 
(.J = is symmetric about this axis. 



as the solution for 7 



CO 




z [Mpc] 



Figure 3. Density deviation profile along the z-axis at difi'er- 
ent times for an initial density perturbation & oc sin(A:z) for 
k = 27r/8Mpc~^ along the z-axis, for 7 = a; = 0. The solid line 
shows the profile of the inhomogeneities today and the dashed line 
shows the same inhomogeneities at a redshift z^^^^ = 5. This 
shows how an initial sinusoidal deviation in the density changes 
its profile at late times due to the non-linear growth of the inho- 
mogeneities. 



grangian 






)\s^{ 



+ 



(16) 

where A is an affine parameter. From this we find the first 
null geodesic equation to be 



SS |(t)' + 



dj/N 



' + z 



Mtr] 



+ZP+US' (t) = 



' dX'-' ' 



(17) 



where an over-dot denotes differentiation with respect to t. 
The second and third equations take the similar form 
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dx 



+ 2^--2BZp+{x + ^) - 



dA2 



S dX 



dz 



and 



d y , r,-dx dy 



dA2 



+ 2 



dA 

5 dA 



2BZl3+{y + Lu) 



dX 



= 0, 



(18) 



(19) 



the last null geodesic equation, governing the photon's mo- 
tion along the z-axis, is given by 



a z 
dA2 






2 -((/3+).{/+ + B[(^ + 7)' + (y + ^)1} 
+ 2B/3+[(a; + 7)7, + (y + cj)a;.]) 



2 — 
dA 



dS 
d\ 

S 



^+2J3(x + 7)t + 2i?(i/ + a;)§ 



(20) 



where a subscript z denotes differentiation with respect to 
the Cartesian coordinate z. These differential equations con- 
tain derivatives with respect to three different variables, t, A 
and z. Since only the Iree lunctions are functions of z and we 
will specify them on a case by case basis, we can consider 
these derivatives to be known. We would like to consider 
functions of only one variables and not two, A and t. Since 
we are considering null geodesies here, we can specify A, such 
that 



dX dt' 



(21) 



where E 



u'^la is the energy of the photon. Using this 

"* "^ ^OXt we can 



relationship and a new time variable r — ^Hq 
simplify the set of differential equations to 

E' 



E 



= S" 



{ 



^HiQA 



Z^ + f'+P^ 



I '-' / '2 I '2 

+ ^i^ +y 



)}, 



+ 2 



.5' 



E 



y + 



+ 



+ p+ 



2^ + ^ 
S E 



2BZp+{x + -i)z''- 



y'-2BZp+{y + uj)/- 



(22) 



(23) 



(24) 



2z' 



S 



2 E 



f'+ + 2B{x + -i)x' + 2B{y + uj)y' 



0, (25) 



where a dash denotes differentiation with respect to t. Initial 
conditions here should be chosen according to the situation 
that is to be modelled. We will always integrate starting 
from the observers position (which we denote by O), which 
we will therefore place at the origin, x\o = y\o ~ z\o ~ 
as the position, where the observer is situated, E\o = 1 
as this is just a normalisation, i.e. the redshift in an inho- 
mogeneous universe z^^ = E/E\o — 1 and thus E\o ~ i 
means that z''^ = E — I. The initial conditions x'\o, y'\o 
and z'\o are chosen depending on into which spatial direc- 
tion we would like to perform the light tracing. These null 



geodesic equations are the most general that we can derive 
in the given space-time. However, our metric only allows 
us to freely choose the initial matter distribution along the 
z-axis (for 7 = a; = 0, otherwise along a path dictated by 
those two functions). Therefore, we will mostly be interested 
in the propagation of light rays along this 'special' z-axis. 
Hence it would be of interest to investigate to what degree 
the above differential equations simplify, if we only consider 
the case where 7 = cj = and light rays only travel along 
the 2-axis. In this case we find the much reduced system of 
differential equations 



£1 
E 



S 



F' 



and 



(26) 



(27) 



where we have used the null constraint from the line element 
to reduce the order of (|25p . We can use these differential 
equations to trace single photons into the past, finding their 
position and energy at any given cosmic time t. However, to 
be able to plot the Hubble diagram, we also need information 
about how bundles of light rays behave in this space-time 
and hence we need to consider the Sachs optical equations. 



4 THE SACHS OPTICAL EQUATIONS 

To describe the evolution of a bundle of light rays one needs 
to specify its expansion 6, shear a, and rotation oj, which 
are the quantities whose evolut ion is described by the Sachs 
optical equations (|Sachslll96ll ). H owever, in this analysis , 
we closely follow the notation of IChandrasekhaiJ ([1993), 
who put the optical scalar equation s in the context of th e 
Newman-Penrose formalism, see also lStephani et al.l l|2003r ). 
Since we are considering point-like sources, we can ignore the 
rotation uj of the light bundles. Hence, the Sachs equations 
take the form 

^ + e' + \af = ^oo 

da „ -, 

— + 2^6) = *o, 
dX 

where 



and 



— Rabl I 



*o = -Cabcdl"''rn'l''m'' 



(28) 
(29) 

(30) 

(31) 



Here ^Pq and 0oo are, respectively, the zeroth Weyl scalar and 
one of the Ricci scalars of the Newman-Penrose formalism 
and represent the Weyl focusing and Ricci focusing in the 
direction of Z" . Rat is the Ricci tensor, Cabcd the Weyl tensor, 
Z" is the affinely parametrised tangent vector to the null 
geodesic defined as 



dx°_ 
dX 



(32) 



and m" is a complex vector that is orthogonal to Z", null 
and has magnitude of 1. The two vectors I"' and m" are part 
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of a complex Newman-Penrose canonical null tetrad. The 
expansion 9 and shear o are precisely 






ka:b)^ 



(33) 

(34) 



We emphasise that for the moment Z" points in a generic 
direction a nd therefore is not the same Z" as in our previ- 
ous paper (JMeures fc Brunill2011^ : however, the advantage 
of using the canonical Newman-Penrose formalism, as pre- 
sented in Chandrasekhar (1992), is that it allows us to easily 
express *o and 4>oo in the equations above in terms of Weyl 
and Ricci scalars in the special null tetrad adapted to our 
metric (see below). The expansion of the bundle of light rays 
is not a direct observable though and so we would rather like 
to consider the angular diameter distance dA and the lumi- 
nosity distance di. One finds the relation 



d\ 
dA 



(35) 



and Etherington's theorem (|Etheringtonlll933l ) states that 

dL = {l + z"'fdA, (36) 



where z is the redshift in a general, inhomogeneous uni- 
verse. In terms of the angular diameter distance, the Sachs 
equations take the form 



dX^ 
da 
d\ 



+ 



Poo — Fj 

d{dA) 
, d\ ^ . 



dA 



= *o. 



As in Sec. [21 we will use the time variable r = |j/oVf^Ai 
and we also introduce 



(39) 



^Jml^K 



and hence we can write the Sachs equations in the form 



C^A +d!. 



E 



<^oo 



9 E^H^Qa 3 E^ 

2 
dA" SV^EH^ 



dA, 



Z——a 



*o 



(40) 
(41) 



where, again, a dash denotes difi'erentiation with respect to 
the time variable r. To calculate the Ricci focusing term 
000 and the Weyl focusing term ^o, we need the form of the 
complex null tetrad. We shall firstly consider the case where 
the photon travels along the z-axis and then generalise the 
result to light rays travelling in any direction. For a light ray 
travelling along the z-axis we find 



and 



^"=£(1,0,0,—) 



7!(°'i'i'°)' 



(42) 



(43) 



where E, again, is the energy of the p hoton. This tetrad is 
very similar to the one derived in i Meures fc Bruni (2011 ) , 
except for the E factor in Z" and so, using this tetrad, we 



find the only non-zero Weyl scalar to hqj ^2. Hence, for light 
bundles along the z-axis, we find the two focusing terms 

1 



000 



:E''p, 



and 



while 



*o 



*2 = -pS 
6 



(44) 



(45) 



(46) 



This brings the Sachs optical equations along the z-axis into 
the form 



dA + dA-rr = 



9 Hi^A 



3 £2 



dA, 



+ 



2^<J = 0. 
dA 



(47) 



(48) 



This system of equations should be integrated from today 
back into the past and so we need to set initial condi- 
tions today, say to and we set a\o ~ 0, dA\o ~ and 
d^lo = —2/{3E\oHo\/Qa)- Given these initial conditions, 
it is apparent from Eq. (|48p that a\o ~ implies the trivial 
solution cr = and hence we only have to consider Eq. (|47|l 
and the initial conditions associated with it. Essentially this 
means that along the z-axis the Weyl focusing is zero and 
hence the light bundles do not experience any shear, how- 
ever, the Weyl focusing is an effect we are interested in, as 
it might have non-negligible effects on the angular diameter 
distance. Hence, we will now generalise the above treatment 
to be able to consider light bundles that do not travel along 
the z-axis. 

We are interested in the case where j — u — 0, since 
these two functions only displace the center of deviations, 
and therefore in this special case, o ur model displays an ax - 
ial symmetry about the z-axis, see lMeures fc Brunil ||201H ). 
This implies that considering light rays in the j/-z-plane is 
completely general, as one could always do a rotation about 
the z-axis without changing the metric but making the tan- 
gent vector point out of the j/-z-plane. We name the angle 
that the tangent vector Z" subtends with the z-axis a and 
all quantities in the rotated system are denoted with a tilde. 
The null tetrad in the rotated frame then takes the form 



l" = £(1,0, 



sin(Q) cos(q) , 



S 



sz 



and 



V2 



(0, 



1 



-zcos(a) — isin(a). 



SZ 



(49) 



(50) 



Since we are dealing with pure ACDM, dust and a cosmolog- 
ical constant, it follows from EFEs that the Ricci focusing 
term does not depend on rotations in the basis vectors and 
so we find that 



— E p. 

2 ^ 



(51) 



^ In doing this, we have chosen a null tetrad which is especially 
adapted to the metric: having 'I'2 as the only non-zero Weyl 
sca lar is characteris t ic of the Petrov type D of our space-time, 
see lMeures &: Brunil 1 I2OIM for more details. The fact that a met- 
ric o£_the_Jorm_we_are_using_^of Pertov type D was first shown 
bv lBarnes fc RowlinesonI l ll989l) . 



© 2011 RAS, MNRAS OOO.nHTSl 



Redshift and distances in inhomogeneous ACDM 7 



However, deriving the Weyl focusing term in the rotated 
frame is not as straight forward and we need to consider the 
effect of rotations in the complex null tetrad on the Weyl 
scalars. In the complex null tetrad constructed for light rays 
travelling along the z-axis, we found that the only non-zero 
Weyl scalar was ^2, Eq. (|46p . Rotating the original com- 
plex null tetrad to coincide with the the physical situation 
of light propagation at an angle a with the z-axis requires 
four separate canonical rotations in the complex null tetrad, 
see Appendixl^for the details. These rotations have the ef- 
fect of making all five Weyl scalars non-zero in general, but 
expressible in terms of the original ^2, in particular, we find 

*o = -3sin^(a)£^*2 = --sin {a)E'^p5, (52) 

where p is the background density and 5 is the density devi- 
ation. Therefore, we find the general Sachs optical equations 
for our space-time 



di 4- d. 



EL 

E 



9 R^^K 



dA, 



2 — CT = 

dA 



3 £2 



sm 



3^/3 Hl^K 



pS, 



(53) 
(54) 



where we choose initial conditions again as ajo = 0, dAlo = 
and d'j^lo = —2/{3E\oHq\/TIa) and integrate into the 
past. Clearly, choosing 5-\o =0 here does not imply the 
trivial solution for a, since the Weyl focusing term is non- 
zero in general. 



5 RESULTS OF THE LIGHT TRACING 

In this section, we would like to present how the results of 
the light tracing we performed in our model differ from the 
standard FLRW results commonly used. To find the position 
and redshift of the bundle of light rays, we need to integrate 
Eqs. (|22|l - (|25| l in general, whereas, to find the angular diam- 
eter distance and shear of the bundle, we need to integrate 
Eqs. (|53|l and ((54)). Inspecting the last two equations, we 
find that they are coupled to the geodesic equations and 
hence, we need to solve all six differential equations simul- 
taneously. From the angular diameter distance, we can find 
the luminosity distance, using Eq. (|36|) . but we would also 
like to compare the distance modulus we find to the standard 
FLRW result and therefore we use 



AdM = dM - df/-^^ = 51ogio ' 



dr 



(55) 



yd^LRW ^ 

where dM is the distance modulus and dL is the luminos- 
ity distance, where a superscript FLRW denotes the same 
quantity in the FLRW background. To compare the angular 
diameter distance and redshift in our model to the standard 
FLRW result, we choose the definitions 



and 



Az- 



AdA 



dA - dj 



^FLRW 



(56) 



(57) 



which should give the reader an intuitive idea of what the 
fractional difference is between the results we derive here 
and the commonly used FLRW values. 



Having derived the geodesic equations and Sachs optical 
equations and having defined quantities to analyse our re- 
sults, we have to consider which physical situations we would 
like to model. For a detailed discussion of which matter dis- 
tributions are possible and a discussion of sing ularities, see 
the relevant sections in lMeures fc Brunil (|201ll ). Here, how- 
ever, we would like to mention the main points that charac- 
terise the matter distributions we can model. In this paper, 
we concentrate on the 7 = a; = case, which reduces the 
freedom of the model to one function, the space distribution 
of the growing mode /3+, which gives us the freedom to set 
the initial matter distribution along the z-axis, whereas we 
do not have any freedom to set the distributions along the x- 
and y-axis. If we only consider under-densities, (/3+ > 0) we 
do not find singularities in any space-time point. However, 
as long as /3+ is negative in some region, which corresponds 
to an over-density in the same region, pancakes can eventu- 
ally form in the model - as also expected from Newtonian 
gravitational collapse. What is important for the analysis 
presented here is that for over-densities which have not yet 
collapsed by today (which is what we are interested in) , we 
can, at any time in the past, find a region around the z-axis 
which is free from singularities and therefore we can perform 
light tracing in those regions. 



5.1 Single mode density distributions 

As we explained in Sec. (2] the linearity of Eq. ((9} allows 
for the validity of a superposition principle for the metric 
function F. In addition, we only consider the growing mode 
whose spatial distribution is encoded in the function l3+{z). 
Therefore, in this section, we first look at harmonic distri- 
butions of matter along the line of sight, i.e. single mode 
sinusoidal distributions. 

The first question that comes to mind is whether to 
consider a distribution of over-densities, under-densities or a 
combination thereof. As a first analysis, we would like to see 
what the effect of either of those three choices is and hence 
in Fig. [4] we present the redshift and the angular diameter 
distance d^ obtained from a model with only over-densities 
(red lines), only under-densities (blue lines) and compen- 
sated density profiles (black lines). The solid lines are for 
light rays along the z-axis and the dashed lines correspond 
to light rays which travel at an angle of 10 and 40 degrees off 
the z-axis, the 40 degrees lines are always the ones further 
away from the respective solid line. No deviation between 
the different angles is visible for the black line as no differ- 
ence from the FLRW curves are visible at this resolution for 
any angle in the compensated case. In Fig. 3] we have cho- 
sen all inhomogeneities to be periodic on a scale of 8 Mpc, 
the initial conditions have been chosen such that at early 



cos( 



2ttz 
SMpc ) 



5undsr oc cos(gjjjfj) - 1 and 



times, (5o„er oc 1 

5comp oc cos(glj^), where Scomp stands for a density pertur- 
bation which is compensated on the above mentioned 8 Mpc 
scale along the z-axis. The amplitude of the over-densities 
in the compensated and only over-density cases correspond 
to (5 ~ 1 today, whereas the under-densities grow to voids of 
5 ~ —0.3 today. The figures show significant deviations in 
the redshift and distance measure, if we only consider over- 
or under-densities along the line of sight. One might get the 
impression here that for the compensated case there are no 
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Figure 4. Redshift z found in our model, as compared to 
the redshift z^^-^^ the same object would have in an FLRW 
model, top panel, and angular diameter distance dji as a func- 
tion of observed redshift z'^ , bottom panel. We are considering 
only over-densities (top, red lines), only under-densities (bottom, 
blue lines) and compensated density distributions (middle, black 
lines). The deviations are periodic on scales of 8 Mpc along the 
^-axis, the solid lines are for on-axis light rays and the dashed 
lines for off-axis rays (corresponding to 10 and 40 degree devia- 
tions). The dashed lines for the compensated case are not visible 
here, as they are not distinguishable form the on-axis case at this 
resolution. 



deviations from the FLRW values, therefore in Fig. [5] we 
have plotted a zoom in on the very small redshift range of 
Fig. |4] for the compensated case only and a periodic devia- 
tions is clearly visible. 

However, when we make actual observations in the Uni- 
verse, we generally assume that we observe along typical 
lines of sight and that the density deviation along this line 
of sight should average to zero, or at least we expect that the 
ensemble average of the density deviation along many lines 
of sight in different directions should average to zero. Hence, 
we would like to consider in the following matter distribu- 
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Figure 5. Redshift found in our model z'^ versus the redshift 
an object at the same distance from us would have in an FLRW 
model z^^^^ . The solid line is the redshift found from the light 
tracing for a compensated density profile, which corresponds to 
(5^1 today, periodic on a scale of 8 Mpc today. We are only 
plotting very small redshifts here, so the deviations found are 
visible. The solid line corresponds to the redshifts we find from 
the light tracing and the dashed line corresponds to the FLRW 
values, plotted for reference. This is a zoom in to the black line 
on the top panel of Fig. |4] 



tions which we can show to average to zero over some char- 
acteristic distance. This is automatically achieved with the 
harmonic sinusoidal distribution we are dealing with which 
clearly implies a zero average over one period of the func- 
tion. Figs. [SJ [7] and [S] show the density deviation 5, redshift 
deviation Az, angular diameter distance deviation Ad^ and 
the distance modulus deviation Adu for density deviations 
on different scales and of different amplitudes. On all plots, 
the solid lines correspond to hght rays which travelled along 
the 2-axis and dashed lines correspond to off-axis light rays. 
The off-axis Hght rays were directed at angles of 5, 10 and 20 
degrees from the z-axis. To distinguish the lines, one can as- 
sume that the ones that deviate from the solid line the most 
are the ones sent at an 20 degree angle and the ones sent at 
5 degrees are hardly distinguishable from the solid lines. We 
are considering here density deviations on different scales, 
ranging from 8 Mpc to 500 Mpc. However, we do not ex- 
pect density deviations to be of the same amplitude across 
all these scales: in general, while the Universe is very in- 
homogeneous on small scales, observations support the idea 
that deviations from homogeneity get smaller and smaller on 
larger scales, wher e homogeneity is reached at some point 
(jSarkar et al.ll2009r ). Generally, we will be considering larger 
amplitude density deviations on smaller scales and smaller 
amplitudes on larger scales. Despite this, interestingly we 
find deviations in the distance measures and redshift to be 
mostly affected by the larger scale density deviations. Con- 
sidering larger amplitude density deviations on larger scale 
would simply amplify this effect but make the density dis- 
tributions less realistic. Note that in Fig. [7] we are using a 
smaller range in redshift than in the other two figures to 
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Figure 6. The eflFect of inhomogencitics with S Ri 0.1 today, 
periodic on a scale of 500 Mpc. The soUd Une corresponds to 
light rays travelling along the z-axis and the dashed lines to light 
rays travelling at 5, 10 and 20 degrees off-axis respectively. 



Figure 7. The effect of inhomogeneities with 5 ss 1 today, pe- 
riodic on a scale of 100 Mpc. The solid line corresponds to light 
rays travelling along the 2-axis and the dashed lines to light rays 
travelling at 5, 10 and 20 degrees off-axis respectively. 



make it easier to tell apart features of the off- and on-axis 
lines. 

Looking at Figs. [6l[7|and[8l we can see deviations from 
the FLRW background values in all the quantities we dis- 
play, however all these deviations seem to be below the 1% 
level, given the conservative assumptions we made. 

We would like to understand what term in the Sachs 
equations, Eqs. (|28|) and (I29|) . is the main cause of these 
deviations and hence we consider the non-FLRW parts of 
the Ricci focusing, the Weyl focusing and the shear result- 



ing from Eqs. (|53p and ([54}. For the Ricci focusing, we will 
introduce the variable 



A<;/) = 



iFLRW 

Poo 



= -Ue'p-E'p), (58) 



where E is the photon energy in the background FLRW 
model. For the Weyl focusing, there is no background con- 
tribution, hence, we just have to consider the Weyl scalar 
^o; for the contribution of the shear, we will consider \a\^ 
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Figure 8. The effect of inhoinogeneities with 5 s:^ 3 today, peri- 
odic on a scale of 1 Mpc. Here we are only showing the results 
from light tracing along the 2-axis for clarity. Two features are 
visible in the lower three panels, the thickness of the lines shows 
the effect of each single oscillation in the density profile, whereas 
the trends of these thick lines shows the overall integrated effect of 
having large amplitude but small scale density inhomogeneities. 



as this is the term present in the Sachs equations, and again, 
a vanishes in the background given our initial conditions. 

In Fig. [9] we compare the above mentioned variables for 
a representative matter distribution: we choose & periodic 
on scales of 100 Mpc along the z-axis with an amplitude of 
5 ~ 1 today. We have performed the integration for light rays 
which travel at an angle of 20 degrees with the z-axis. From 



"^ 



O 





Figure 9. Plot of the non-FLRW contributions in the Sachs equa- 
tions as a function of redshift. Here we have used compensated 
density deviations on scales of 100 Mpc which would have an am- 
plitude of (5 ft! 1 today and the light tracing was performed at an 
20 degree angle with the z-axis. Where A0 = 0oo — </'od^^^> ^0 
is the zeroth Weyl scalar and cr is the shear of the light bundle. 



the plots we can see that the Ricci focusing deviation dom- 
inates over the Weyl focusing. Although from Eq. (|58|l we 
see that A(^ « Ep5 so that the source of both the Ricci and 
Weyl focusing is p5, the Weyl focusing is strongly suppressed 
by the sin (a) factor in Eq. (|52p . Finally, we find that jaj is 
much smaller than the two focusing terms, so that it gives a 
negligible contribution in Eq. (I28p and, consequently, to the 
deviations from FLRW in the distance measures. 



5.2 Mode coupling and its effects 

The above results show that single mode density deviations 
do not have a large effect on the redshift, angular diameter 
distance and distance modulus for the compensated profiles 
we considered. From this we cannot conclude though that a 
density profile, where the metric function F is the sum of 
many modes, has a small effect as well, since the structures 
in our model grow non-linear ly. In other words, an initial 
deviation consisting of the superposition of two modes in 
F may excite many different modes in 5 during its non- 
linear growth and the resultant density profile might have 
completely different effects on the redshift and the distance 
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modulus. The real Universe clearly does not consist of only 
one wavelength deviation and therefore the step of including 
several modes should make our analysis more realistic. 

In general in our model we find that the growth of sin- 
gle mode density deviations does not depend on their wave- 
length, that is, if we only had deviations periodic on a 1 Mpc 
scale, they would grow to the same amplitude as if we only 
had deviations on a scale of 100 Mpc. If we had an initial 
small density deviation which is a superposition of those two 
modes in F, the growth of 5 would be quite different to the 
individual modes. In Fig.[TO]we show how modes on scales of 
1 Mpc, 20 Mpc and 100 Mpc interact. From these plots one 
can see how peaks in the long wavelength perturbation cause 
the short wavelength perturbations to grow non-linearly and 
it becomes clear that the profiles are not superposition of the 
individual modes any more. More precisely, in the top panel 
of Fig. [To] we show the growth of a 1 Mpc mode together 
with a 20 Mpc mode and their interaction. It is clear that 
the growth of the short wavelength mode gets an extra non- 
linear kick from growing on top of the larger scale mode. On 
the other hand, peaks of the short scale mode that grow in 
the voids of the larger scale mode are depressed. The same 
qualitative behaviour can be observed in the middle panel 
of Fig. IIUI where we now consider the 1 Mpc mode together 
with a 100 Mpc mode, plotting on the same length scale as 
the top panel. In this case the peak of the larger scale mode 
(in the middle of the figure) is very broad and the short 
wavelength mode is growing almost as if on top of a differ- 
ent background. However, non-linearity is again important 
and the peaks are much higher than they would be if sim- 
ply raised by this "new background" . Finally, in the bottom 
panel we show the effect of adding the larger 100 Mpc mode 
to the deviations in the top panel. The red profile for the 
coupled 1 Mpc and 20 Mpc modes of the top panel is shown 
again in red in the bottom panel. The green profile shows 
the effect of coupling the three modes together. This bottom 
plot in Fig.[lO]therefore shows that even adding a small am- 
plitude {5 ~ 0.3 today) 100 Mpc mode to the 1 Mpc and 
20 Mpc allows the peaks on the shorter scale to grow much 
stronger, more than a factor of 4 than without. 

To investigate the effect of the mode coupling on the 
redshift and distances, in Fig. [TT]we show how two coupled 
modes, one 1 Mpc mode and one 100 Mpc mode, affect the 
redshift and distances. The two modes have been chosen to 
be of equal initial amplitude and combine to result in struc- 
ture of (5 ~ 1 today. The amplitude of deviations in redshift, 
angular diameter distance and distance modulus are com- 
pletely dominated by the effects of the long wavelength, 100 
Mpc deviation. This analysis has been done for many more 
pairs of modes, from 1 Mpc to 500 Mpc, and the results al- 
ways seem to be dominated by the long wavelength modes. 
The fact that the main effects in redshift and therefore in 
distances is dominated by the larger scales inhomogeneities 
should not come as a surprise. The basic mechanism at work 
here is the same as in the Rees-Sciama and the integrated 
Sachs- Wolfe effect. In an expanding universe photons travel 
through dynamical inhomogeneities. When the characteris- 
tic scale of the inhomogeneity is negligible compared to the 
Hubble radius, the effect of the expansion is negligible and 
therefore a photon will come out of the inhomogeneity with 
the same energy that it had when it entered it. Instead, in 
going through a large scale inhomogeneity, photons have to 




z [Mpc] 



Figure 10. Illustration of how different modes interact in the 
density profile. Each profile is plotted at a redshift of z^^^^ = 5 
(smaller amplitude curves) and today. The top panel demon- 
strates the interaction of a mode that is periodic on 1 Mpc, blue 
curves, and a mode periodic on 20 Mpc, black curves, and how 
the two modes interact, red curves, if superimposed as initial con- 
ditions. The same is shown in the middle panel for modes periodic 
on 1 Mpc, blue curves, and 100 Mpc, black curves, and their in- 
teraction, red curves. In the bottom panel, we show how an initial 
superposition of the 1 Mpc and 20 Mpc modes, red curves, behave 
compared to an initial superposition of all three, 1 Mpc, 20 Mpc 
and 100 Mpc modes, green curves. 



go through a different potential well when they come out of 
the inhomogeneity than when they were entering it, chang- 
ing their energy in the process. 



5.3 Peaks and voids of arbitrary profile 

After having analysed a variety of different sinusoidal sin- 
gle mode matter distributions and their coupling, we would 
like to investigate whether a more complex matter distribu- 
tion would give a more significant deviation from the FLRW 
background. As shown in Fig.|4j choosing a matter distribu- 
tion which is not compensated, i.e. where integrating F over 
any distance along the 2-axis does not give zero, has signifi- 
cant effects on redshift and angular diameter distance. Now 
we would like to investigate whether compensated profiles 
which are more complex than simple single mode sinusoidal 
can have significant effects as well. To this end, we choose 
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Figure 11. The effect of infiomogeneities with 5 fa 1 today, with 
the initial condition being a superposition of two modes, one pe- 
riod on scales of 1 Mpc and the other on 100 Mpc. The initial 
amplitude of the two different perturbations was chosen to be 
the same. Here we are only showing the results from light tracing 
along the z-axis for clarity. 



Figure 12. The effect of inhomogeneities with 5 si 2.3 today, 
witli the choice of the initial profile shown in Eq. II59I I. Here we 
are only showing the results from light tracing along the 2;-axis 
for clarity. This graph clearly shows that choosing non-sinusoidal 
matter deviations does not change the results significantly. 



the initial profile for each over-density to take the form 



5 oc cosh ( 



lOMpc 



)-c, 



(59) 



where C is a constant, and propagate hght rays through a 
periodic array of such shapes. To ensure that this distribu- 
tion is compensated, the constant C needs chosen carefuUy. 
In Fig. [12] the results of this analysis are shown, where we 



have chosen the distance between the peaks to be 100 Mpc 
today. The deviations from the FLRW results are not signif- 
icantly different from the ones obtained by using sinusoidal 
distributions; in particular the dominant parameter in the 
effects on redshift and distances is the maximum length scale 
of the deviations. 
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6 CONCLUSIONS 

In this paper we have analysed the effects of non-hnear 
structure on redshift and distance measure s using the exact 
solut ion developed in our previous paper (jMeures fc Brunil 
I2OIII ). This model, described by the line element in syn- 
chronous comoving coordinates, allows us to choose an ar- 
bitrary matter distribution along one line of sight with the 
growth rate of structure and the density distribution away 
from this a:xis being set by EFEs. A remarkable feature of 
our model is that the inhomogeneities are described by a 
single metric function F, which satisfies the same linear sec- 
ond order differential equation satisfied by the linear density 
perturbation S in Newtonian perturbation theory as well as 
in relativistic perturbations when the synchronous comoving 
gauge is used. Therefore F satisfies a superposition princi- 
ple and, in particular, extends into the non-linear regime the 
same growing and decaying modes that S shows in the linear 
regime. 

We have developed the null geodesic equations and the 
Sachs optical equations in our model for light rays travel- 
ling in arbitrary directions. This set-up has then been ap- 
plied to different physical situations, considering single si- 
nusoidal mode deviations in the density, the coupling of two 
and three harmonic modes as well as a more complex matter 
distribution described by an array of peaks and voids along 
the line of sight. Furthermore, we have investigated which 
terms in the Sachs optical equations are mainly responsi- 
ble for deviations from the FLRW values. Additionally, we 
have analysed the interaction of two and three modes in the 
growth of structure within our exact non-linear framework. 

We consider the redshift and distance measures for sin- 
gle mode density deviations on different length scales and of 
different amplitudes. The largest effect for the redshift and 
distance measure to deviate from the FLRW results seems to 
be obtained for larger density deviations and larger scales. 
The results are shown in Figs. [6l [7] and [S] Given our con- 
servative assumptions on the density profiles, all deviations 
are below the 1% level. 

Even if the metric function F satisfies a linear equation, 
our model is non-linear and so the effects of two modes on 
the redshift and distances does not need to be the same as 
their combined effects, hence we considered how modes in- 
teract in the growth of the density deviations and how the 
redshift and distance measures are affected by them. We 
find that combining small and large scale density deviations 
has significant effects on the growth of structure, in that 
the peaks of the large scale modes significantly enhance the 
growth of small scale deviations, see Fig. [10] for plots illus- 
trating this point and the text in Sec. 15.21 For more details 
on the growth of structure in our model, see lMeures fc Brunil 
(|201ll ) and for a more general discussion of the relativistic 
behaviour of inhomogeneities a nd overall proper t ies of the 
class of solut ions we consider, see lMatarrese et al.l (|l994al lbr): 
iBruni et al.l (ll995a,b). The interaction between short and 
long wavelength modes observed in our model does not seem 
to have a significant effect on the redshift and distances mea- 
sures though, as the effect of the long wavelength deviation 
remains dominant despite the presence of small scale devi- 
ations, see Fig. [TT]for the results of a 1 Mpc and 100 Mpc 
(today) mode combination. This implies that mode coupling 



does not provide significantly larger deviations on the red- 
shift and distances than the individual modes. 

To generalise our result to density distributions where 
peaks and voids are more pronounced than in a single mode 
sinusoidal, we have considered an array of density profiles 
given in Eq. (|59|) . which provides quite peaked over-densities 
and large voids separating them, choosing a typical array 
scale of 100 Mpc today. The results of the light tracing for 
this distribution is given in Fig. 1121 where the density distri- 
bution is shown in the top panel. Given this density profile, 
which is clearly non-sinusoidal for all times, the deviations 
from the FLRW redshift and distance measure are still com- 
parable in amplitude to the results found in the single 100 
Mpc mode analysis, with the large scale deviations being 
below the 1% level. 

For all these different density distributions, we inves- 
tigated which of the terms in the Sachs optical equations 
is dominant in providing the deviations, see Fig. [5] for the 
different terms for a single mode deviation. We find that in 
all cases the Ricci focusing term is dominating, while the 
effect of the shear on the angular diameter distance seems 
to be vanishingly small. Due to its special geometric char- 
acter (the space-time is Petrov type D), in our model the 
Weyl focusing is exactly zero along the z-axis of symmetry 
and it is sub-dominant with respect to the Ricci focusing in 
direction at an angle a with respect to the z-axis because of 
a sin^(a) factor. However, both the Ricci and Weyl focusing 
are proportional to the density deviation S and so we may 
expect that in a more general space-time they would be of 
the same order. 

All the above mentioned results are for density devia- 
tions which we refer to as compensated, i.e. where the met- 
ric function F averages to zero along the line of sight at 
all times, so that the initially small density profile is also 
compensated. The results from considering density profiles 
which do not average to zero initially, however, are quite dif- 
ferent. In Fig. |4j we show that over-dense and under-dense 
lines of sight have significant effects on the redshift and an- 
gular diameter distance. This implies that if we, on average, 
observe along lines of sight which are more or less dense than 
the background, we may need to expect significant effects. 
Therefore it is important to understand whether the lines of 
sight we observe along are really average "skewers" through 
the Universe matter distribution. It also emphasises the im- 
portance of identifying the correct background density with 
observations, as an over- or under-estimate may affect our 
interpretation of observational results significantly. Under- 
standing these effects is crucial to our interpretation of cos- 
mological observations and hence we leave a deeper analysis 
for future work. 

Finally, part of the motivation for our work has come 
fr om the strong devia t ions fro m the standard FLRW results 
in IClifton fc Ferreiral (|2009al ) , and therefore we would like 
to briefly compare our results here. Considering density dis- 
tributions which are initially compensated along the line of 
sight cannot provide deviations in the redshift and distance 
measures as large as found in IClifton fc Ferreiral (|2009al ). 
This does seem to be in agreement with most Swiss-Cheese 
type analyses. However, considering lines of sight which have 
an average density which is lower than the bac kground den- 
sity can provide similar results as found i mClifton fc Ferreiral 
(|2009ah . 
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APPENDIX A: TETRAD TRANSFORMATIONS 

Here we would like to present how to find the Weyl focus- 
ing term ^o for light rays travelling at an angle a from 
the z-sxis. In the Newman- Penrose formalism, the five Weyl 
scalars ('I'o-'I'4) are contractions of the Weyl tensor with a 
complex null tetrad, Z °, n", m" and m'^; for our notation see 
Chandrasekhan (|l992h . In our first paper, (IMeures fc Brunil 
201ll 'l. we have derived the Weyl scalars for the null tetrad 



r = 






(Al) 
(A2) 
(A3) 



where S and Z are the metric functions and m" is simply 
the complex conjugate of m°. Please note that we are using 
slightly different notation he re than in the fir s t pap er, as we 
follow the exact notation of IChandrasekharl ([1993) here to 
avoid confusion when referring to this book. In this special 
null tetrad, the only non-zero Weyl scalar is ^2. In this 
paper, we are interested in light tracing and we find that the 
tetrad vector /° = da;°/dA and so depending on the direction 
the light rays travel in, the vector /" is going to change. 
Hence we need to understand how changes in the complex 
null tetrad affect the Weyl scalars. Given the properties of 
complex null tetrads, there are only three distinct types of 
transformations, those of type I 



I ^ I , m 
and n" 

of type II 

n ^ n , m - 
and r 

and of type HI 

r^A-^r, n" 

and m" 



771 + al ,ni — )-m + a I , 
n + a m + am + aa l , 



(A4) 



m + on ,m. -^ m + n , 

r +h''ni' +hrn' + bh'n'', (A5) 



An ,m -^ e m 



— iO — a 

e m 



(A6) 



where a and b are complex and A and 6 are real valued 
functions. Each of these three rotations h as the effect of 
mixin g the Weyl scalars in a certain way, see 1 Chandrasek han 
(|l992r ) for the exact relations. To find the tetrad determined 
by Eqs. (|49p and H50[) from the above tetrad, we need to 
perform a combination of these transformations. We have 
used, in the given order, a transformation of type II with 



b = i- 



3(a) 



sin(Q:) 

a transformation of type HI with 

» ,,1 ^ cos(q) , ^ 

Al = 2 „// and 9^ = 0, 

sm^(a) 



(A7) 



(A8) 



a transformation of type I with 

.cos(a) 



1 



sin(a) 
and finally a transformation of type III with 

^ and 62 = 0. 



^/2E 



(A9) 



(AlO) 



From these rotations, we obtain the null tetrad given in Eqs. 
(gg and dnSJ from the null tetrad in Eqs. fSTjl - (|M)) . Given 
the transformation rules of the Weyl scalars for the above 
rotations, we find all five Weyl scalars to be non-zero in 
general and specifically, we find that 



*o 



-3sin^(a)£'^*2, 



(All) 



where the tilde denotes the quantity after the rotations. This 
is the Weyl focusing term for light rays travelling at an angle 
a to the 2-axis. 
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